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DEFENCE  RESEARCH  ESTABLISHMENT  SUFFIELD 
RALSTON,  ALBERTA 

SUFFIELD  TECHNICAL  NOTE  NO.  444 

DESIGN  OF  A  FILM-COOLED  ENTRAINING  DIFFUSER  (U) 


by 

S.B.  Murray 


ABSTRACT 

A  film-cooled  entraining  diffuser  is  described  which  consists 
of  a  series  of  staged  cylindrical  rings,  each  overlapping  the  adjacent  one 
so  as  to  create  annular  slots  for  the  entrainment  of  surrounding  ambient  air. 

A  two-step  iterative  design  procedure  is  outlined.  In  step  one  an 
analysis  similar  to  that  first  employed  by  von  Karman  is  used  to  calculate  the 
rate  at  which  air  is  drawn  into  each  slot.  In  step  two  these  flow  rates  are 
used  in  a  downstream-marching,  iterative,  implicit  finite-difference  method  to 
calculate  the  development  of  wall  jet  boundary  layers  downstream  of  the  slots. 

Details  about  the  design  and  manufacturing  of  a  three-ring  model 
structure  are  presented  and  proposed  future  experimental  validation  is  dis¬ 
cussed. 

(U) 
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is  a  cross-sectional  area  of  flow, 
is  a  specific  heat  at  constant  pressure. 

is  the  diameter  of  the  duct  to  which  the  film-cooled 
entraining  diffuser  is  fastened. 

is  the  diameter  of  the  n-th  ring  element  at  the 
downstream  end  of  the  film-cooled  entraining  diffuser. 

is  a  radiation  shape  factor. 

is  the  number  of  an  interior  nozzle  ring  element. 

is  a  loss  coefficient  associated  with  the  entrance 
to  an  entrainment  slot. 

is  the  length  of  a  nozzle  ring  element,  or  i  is  a 
mixing  length  in  the  turbulence  model. 

is  the  overall  length  of  the  film-cooled  entraining 
diffuser. 

is  the  number  of  nozzle  ring  elements  comprising  the 
film-cooled  entraining  diffuser. 

is  the  static  pressure  in  the  boundary- layer 
calculations  and  is  assumed  to  be  a  function  of  x  alone. 

is  the  uniform  static  pressure  at  the  entrance  or  exit 
of  a  nozzle  ring  element. 

is  the  molecular  Prandtl  number. 

is  the  turbulent  Prandtl  number. 

is  the  local  radius  in  axisymmetric  flow. 

is  the  gas  constant  for  primary  and  entrained  streams. 

is  the  fluid  static  temperature. 

is  the  streamwise  (x-direction)  component  of  velocity, 
is  the  transverse  (y-direction)  component  of  velocity, 
is  a  slot  width. 

is  the  streamwise  coordinate,  measured  parallel  to  the 
wall . 

is  the  transverse  coordinate,  measured  normal  to  the  wall, 
is  the  fluid  eddy  thermal  conductivity. 

is  the  boundary- layer  thickness. 

is  von  Kantian' s  mixing-length  constant. 

is  a  constant  relating  mixing  length  to  boundary-layer 
thickness  in  the  outer  region  of  the  boundary  layer. 


Nomenclature  (continued) 


u  is  the  fluid  dynamic  viscosity. 

vt  is  the  fluid  eddy  viscosity. 

II  is  Coles'  law  of  the  wake  profile  parameter. 

p  is  the  fluid  static  density. 

Subscripts 

e  denotes  value  in  the  free  stream. 

1  denotes  value  for  an  arbitrary  nozzle  ring  element. 

w  denotes  value  at  the  wall. 

!  denotes  value  in  the  main  stream  at  the  entrance 

to  a  nozzle  ring  element. 

2  denotes  value  in  the  entrained  stream  at  the 
entrance  to  a  nozzle  ring  element. 

3  denotes  value  in  the  combined  stream  at  the  exit 
of  a  nozzle  ring  element. 

Superscripts 

k  is  zero  for  two-dimensional  flow  and  unity  for 

axi symmetric  flow. 

'  denotes  a  time  fluctuating  quantity. 
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1.0  INTRODUCTION 


For  over  a  decade  now  common  use  has  been  made  of  thrust  augmenters 
to  increase  the  propelling  force  of  aircraft  gas  turbine  engines,  particularly 
at  low  speeds  or  in  a  static  situation.  In  its  simplest  configuration  a 
thrust  augmenter  takes  the  form  of  a  constant  diameter  cylinder  or  nozzle 
which  is  positioned  concentric  to  and  slightly  overlapping  the  exhaust  duct  as 
shown  in  Figure  1.  The  action  of  high  viscous  shear  forces  and  turbulent 
entrainment  act  to  draw  surrounding  ambient  air  through  the  overlap  region  and 
into  the  nozzle  where  it  mixes  in  a  turbulent  fashion  with  the  engine  exhaust 
gas.  Although  the  present  study  takes  advantage  of  precisely  this  phenomenon, 
it  does  not  concern  itself  with  thrust  augmentation,  but  rather  with  free 
entrainment  of  ambient  air  for  the  purpose  of  achieving  a  physically  compact 
diffuser,  efficient  exhaust  stream  dispersion,  and  efficient  wall  film  cooling, 
or  any  combination  of  these  possibilities. 

In  an  industrial  application,  for  example,  where  maximum  output  of 
shaft  power  may  be  of  most  concern,  it  is  beneficial  to  realize  a  pressure 
recovery  by  diffusing  the  turbine  exhaust  stream  before  discharging  it  to  the 
atmosphere.  Such  a  process  could  be  accomplished  by  the  use  of  a  simple 
diffuser,  but  this  is  often  impractical  due  to  the  small  angle  of  divergence 
imposed  on  the  constraining  walls  in  order  to  prevent  flow  separation.  A  much 
higher  effective  rate  of  diffusion  is  possible  with  a  slotted  diffuser  as 
shown  experimentally  by  Frankfurt  (1975)  and  others.  This  is  so  because  the 
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FIGURE  1  :  A  TURBOJET  ENGINE  AND  SIMPLE  THRUST 

AUGMENTER, 
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injected  air  streams  form  a  series  of  energetic  wall  jet  boundary  layers  which 
act  to  delay  the  onset  of  flow  separation  in  a  manner  similar  to  that  of  slot 
blowing  over  the  upper  surface  of  airfoils. 

In  another  conceivable  application  rapid  cooling  and  dispersion  of 
the  exhaust  plume  may  be  necessary  for  reasons  of  safety  or  military  strategy 
or  to  minimize  the  effects  of  local  thermal  pollution.  In  such  a  case  one  of 
the  design  objectives  would  be  to  incorporate  relatively  large  slots  in  order 
to  promote  entrainment  of  the  substantial  volumes  of  ambient  air  needed  to 
lower  the  bulk  temperature  of  the  combined  stream  at  exit. 

The  objective  in  yet  another  application  may  be  to  isolate  the  dif¬ 
fuser  walls  from  corrosive  or  hot  exhaust  gas  either  for  reasons  of  military 
strategy  or  to  prevent  structural  damage  resulting  from  differential  thermal 
expansion  or  surface  pitting. 

In  any  event  it  is  evident  that  criteria  governing  the  design  of  a 
film-cooled  entraining  diffuser  will  be  dictated  by  its  end  use. 

2.0  THE  DESIGN  PROCEDURE 

Given  the  objectives  of  a  particular  diffuser  its  design  must  be 
carried  out  subject  to  certain  geometric  constraints.  For  example,  it  must  be 
compatible  with  an  exhaust  duct  of  diameter  Dj  and  maintain  an  overall  length 
of  L  and  an  exit  diameter  of  D2. 

Inputs  to  the  problem  include  parameters  which  characterize  the  gas 
flow  within  the  duct,  such  as  its  mean  velocity  Uj  and  absolute  temperature 
Tj.  If  radiative  heat-transfer  processes  are  to  be  included  in  the  analysis 
the  effective  gray-body  emissivities  of  the  duct  material,  diffuser  material 
and  exhaust  gas  must  be  specified.  Ambient  pressure  and  temperature  are  also 
required  inputs. 

Once  the  constraints  and  inputs  have  been  supplied  the  "design  pro¬ 
cedure"  includes  determination  of  several  more  parameters  including  the  number 
of  nozzle  ring  elements  or  slots  n,  the  length  of  each  ring  element  ,  and 
the  width  of  each  slot  w...  These  parameters  are  illustrated  in  Figure  2. 
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There  are  two  major  steps  involved  in  the  computational  procedure. 
In  Step  1  the  rate  of  entrainment  into  each  of  the  n  slots  is  computed.  Step 
2  utilizes  these  flow  rates  in  order  to  calculate  the  development  of  the  wall 
jet  boundary  layer  downstream  of  successive  slots.  Details  of  each  step  are 
laid  out  below. 

2.1  Step  1:  Calculating  the  Rates  of  Entrainment 

As  mentioned  in  the  introductory  remarks,  it  is  possible  to 
induce  a  secondary  flow  of  ambient  air  by  directing  an  exhaust  jet  into 
a  nozzle  or  series  of  nozzle  ring  elements.  Figure  3  illustrates  how 
secondary  air  is  drawn  into  the  low-pressure  throat  region  of  a  simple  jet- 
nozzle  combination  under  the  influence  of  viscous  shear  forces  and  turbulent 
entrainment.  Here  it  mixes  with  the  exhaust  gas  to  form  a  turbulent  free 
shear  layer.  On  leaving  this  zone  the  flow  undergoes  a  pressure  rise  to 
ambient  static  pressure  as  it  arrives  at  the  exit  plane  of  the  nozzle. 

The  original  calculations  pertaining  to  this  phenomenon  were  per¬ 
formed  by  von  Karman  (1949)  with  a  view  to  increasing  the  thrust  of  turbojet 
engines.  Although  adequate  information  relating  to  thrust  augmentation  is 
available  in  the  open  literature  there  is  an  absence  of  data  with  regard  to 
the  optimization  of  entrainment  systems  from  the  points  of  view  of  plume  and 
wall  cooling.  Consequently,  in  order  to  appreciate  the  potential  of  entrain¬ 
ment  schemes,  the  original  analysis  on  thrust  augmentation  has  been  repeated 
but  with  the  inclusion  of  the  energy  equation  and  with  particular  emphasis  on 
entrainment  performance.  The  importance  of  heat-transfer  effects  is  made 
clear  by  the  theoretical  analysis  of  Quinn  (1976)  who  showed  that  increasing 
the  temperature  of  the  primary  fluid  degrades  the  performance  of  ejectors. 

Briefly,  the  present  analysis  is  as  follows.  Consider  the  simple 
jet-nozzle  combination  of  Figure  3.  The  cross-sectional  area  of  the  nozzle 
is  taken  to  be  constant  and  equal  to  A3  with  pressure  P3  at  the  nozzle  exit 
plane  equal  to  the  undisturbed  ambient  static  pressure.  The  exhaust  jet 
enters  the  nozzle  with  uniform  velocity  Uj  and  temperature  T]  through  area 
Aj.  The  secondary  flow  of  cooling  air  is  through  area  (A3-Aj)  and  with 
uniform  velocity  u2  and  temperature  T2.  It  is  assumed  that  the  exiting  flow 
is  completely  mixed  and  with  monotonic  velocity  and  temperature  profiles  of 
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values  U3  and  T3,  respectively.  Furthermore,  Bernoulli's  equation  can  be 
written  for  the  secondary  flow  between  ambient  infinity  and  a  point  in  the 
throat  just  upstream  of  the  mixing  zone.  Friction  and  turbulence  losses  in 
the  nozzle  are  assumed  to  be  negligible  and  the  nozzle  wall  is  considered  to 
be  adiabatic.  The  governing  equations  for  incompressible  mean  flow  in  this 
system  are: 


Conservation  of  Mass 

p  A  u  +  p  (A  -A  )u  =  p3  A3  u3 

111  2312 

Conservation  of  Momentum 

P  A  -  P  A  =  p  A  u  (u  -u  )  +  p  (A  -A  )u  (u  -u  ) 

23  33  11131  231232 


Conservation  of  Enerc 


P  A  u  (C_  T  +  -i-)  +  p  (A  -A  )u  (C  T  +  -1-) 
lllP.l  ^  2312  P„  2  *■ 


P  A  u  (C  T  +  -5—  ) 


Bernoulli's  Equation  for  the  Entrained  Stream 


P  =  P  + 

3  2 


0) 

(2) 

(3) 

(4) 


Equation  of  State 


P  =  pRT 

where  A  is  cross-sectional  area  of  flow, 

u  is  streamwise  velocity, 

T  is  static  temperature, 

P  is  static  pressure, 

p  is  static  density, 

K  is  an  entrance  loss  coefficient  associated  with  the 

secondary  stream, 

Cp  is  the  specific  heat  at  constant  pressure,  and 
R  is  the  gas  constant. 


(5) 


The  solution  to  this  system  of  equations  for  the  case  where  the 
primary  gas  is  air  and  where  the  secondary  stream  suffers  no  entrance  loss 
(K=0)  is  shown  graphically  in  Figure  4.  The  ratio  of  jet  area  to  nozzle  area 
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FIGURE  4  :  PERFORMANCE  OF  THE  SIMPLE  JET-NOZZLE 

COMBINATION. 


t 
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is  plotted  along  the  abscissa  while  the  ratio  of  entrained  air  velocity  to  jet 
velocity  is  scaled  along  the  ordinate.  Calculations  have  been  performed  for 
several  values  of  absolute  temperature  ratio,  defined  as  the  ratio  of  exhaust 
gas  temperature  to  ambient  air  temperature.  The  curves  indicate  that  the 
velocity  of  the  entrained  stream  decreases  for  a  given  jet  velocity  as  either 
the  area  ratio  or  the  temperature  ratio  increases.  Particular  attention 
should  be  paid  to  the  gradient  of  the  curves  at  high  values  of  area  ratio. 

The  steep  slopes  imply  that  the  velocity  of  entrained  air  is  very  sensitive  to 
slot  width  in  this  range  and,  since  the  cooling  capability  of  the  secondary 
stream  is  coupled  to  this  velocity,  one  would  anticipate  a  pronounced  variation 
in  wall  cooling  effectiveness  as  the  slot  width  is  varied.  It  should  be 
emphasized  that  although  these  calculations  yield  ideal  results  McCormick 
(1969)  has  shown  that,  for  the  range  of  Reynolds  numbers  and  large  area  ratios 
typical  of  the  present  application,  excellent  agreement  between  calculated  and 
measured  performance  has  been  observed  for  a  variety  of  thrust  augmenters. 

The  unknown  quantities  in  the  system  of  equations  above  include  the 
velocity  of  the  entrained  stream  and  the  temperature  and  velocity  of  the  mixed 
flow  which  exits  the  nozzle.  Although  the  solution  is  relatively  straight¬ 
forward  for  the  simple  jet-nozzle  combination  of  Figure  3,  the  matter  is 
complicated  somewhat  when  several  nozzle  ring  elements  are  staged  together  to 
form  a  multiply-slotted  assembly  similar  to  that  depicted  in  Figure  2.  Under 
these  conditions  the  static  pressure  at  the  exit  plane  of  any  interior  ring 
element  i  is  less  than  atmospheric  pressure  due  to  the  presence  of  the  down¬ 
stream  neighbouring  ring  element  i+1.  The  actual  drop  in  static  pressure, 
however,  is  linked  to  the  rate  of  entrainment  into  slot  i+1,  hence  coupling 
the  performance  of  adjacent  ring  elements.  This  dependency  requires  that  the 
governing  equations  for  each  of  the  n  ring  elements  be  solved  simultaneously, 
an  undertaking  of  some  magnitude  due  to  the  non-linearity  of  the  equations. 

An  alternative  to  this  method  of  solution  and  the  one  employed  here 
is  one  of  iteration.  The  basis  of  such  an  iterative  procedure  is  as  follows. 
The  velocity  of  the  entrained  air  stream  in  the  first  slot  is  guessed  and 
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the  system  of  equations  governing  mean  flow  in  the  first  ring  element  is 
solved  to  yield  the  static  pressure,  static  temperature  and  velocity  of  the 
mixed  stream  at  the  exit  plane.  These  exit  conditions  become  the  entrance 
conditions  for  the  second  ring  element  and,  since  the  static  pressure  at  this 
location  has  just  been  specified,  the  velocity  of  the  entrained  stream  in  the 
second  slot  is  fixed.  This  enables  the  governing  equations  for  mean  flow  in 
the  second  ring  element  to  be  solved  as  they  were  for  ring  number  one.  This 
procedure  is  repeated  for  all  n  ring  elements.  A  solution  is  realized  if  the 
calculated  static  pressure  at  the  exit  plane  of  the  final  ring  element  is 
equal  to  the  undisturbed  ambient  static  pressure,  as  it  must  for  incompressibl 
flow.  A  Newton  root-finding  technique  is  used  to  achieve  convergence. 

A  computer  program  'ENTRAIN'  which  calculates  the  rates  of  entrain¬ 
ment  into  each  slot  of  an  n-slot  assembly  appears  in  Appendix  A.  Up  to  twenty 
ring  elements  can  be  included  in  the  analysis.  The  mean  velocity  and  tempera 
ture  of  the  main  stream,  the  exhaust  duct  and  ring  radii,  the  entrance  loss 
coefficients  and  the  pressure  and  temperature  of  the  reservoir  from  which 
each  slot  draws  its  air  are  the  only  required  inputs.  An  example  will  be 
presented  in  Section  3.0. 

2.2  Step  2:  Calculating  the  Development  of  the 
Wall  Jet  Boundary  Layers 

Once  the  mean  flow  rate  into  each  of  the  n  slots  has  been  deter¬ 
mined  by  the  procedure  in  Step  1  development  of  the  wall  jet  boundary  layers 
in  the  downstream  direction  can  be  readily  computed.  In  the  present  study 
this  is  accomplished  by  numerical  solution  of  the  two-dimensional  or  axisym- 
metric  turbulent  boundary  layer  equations  using  a  downstream-marching, 
iterative,  implicit  finite-difference  method.  Since  full  details  of  the 
model  are  given  elsewhere  by  Murray  (1979)  only  the  highlights  will  be 
presented  here  for  the  sake  of  completeness. 

The  governing  boundary- layer  equations  for  incompressible  turbulent 
flow  in  terms  of  time-averaged  mean  flow  quantities  are: 

Continuity 

|^-  (rkpu)  +  |y  (rkpv)  =  0  (6) 
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x-Momentum 


,  au  ,  3U  _  1  dp  ,  1  3  k#  3U  -t-ttn 

ua3T +  va y "  ‘  p  +  ;r  37  -  pu  v  ) 


Energy 


_ ! _ 1  rk(_JL  H  _ 

p  dv/  r  'Pr  3y  p  v  ' 


State 


p  =  pRT 


(7) 


(8) 


(5) 


with  the  following  Boussinesq  eddy-diffusivity  assumptions  for  the  Reynolds 
stress  and  heat-transfer  terms: 


— 1 — r  3u 

=  pvt  ay 


-pT'v1  -  pa 


3T 
t  3y 


(9) 

(10) 


where  the  definition  of  the  turbulent  Prandtl  number  is 


Prt  -  (11) 

Here  vt  and  at  are  the  eddy  viscosity  and  eddy  thermal  conductivity. 

The  exponent  k  is  equal  to  zero  for  plane  flow  and  equal  to  unity 
for  axisymmetric  flow.  The  coordinate  system  is  curvilinear  in  which  x  and  y 
are  distances  along  and  normal  to  the  body  surface,  with  u  and  v  the  velocity 
components  within  the  boundary  layer  in  the  x-  and  y-  directions,  respect¬ 
ively. 

The  boundary  conditions  associated  with  the  above  equations  for  the 
present  application  are: 

Momentum 


u(x,0)  =  0 
v(x,0)  =  0 
lim  u(x,y)  =  u  (x) 

y-Ko 


(12) 


UNCLASSIFIED 


UNCLASSIFIED 


8. 


Energy 


lim 


y-*® 


3T  /  n\  -  31 

ay  ^X’0)  =  ay  w 
T(x,y)  =  Te(x) 


(13) 


where  subscripts  w  and  e  denote  conditions  at  the  wall  and  at  the  outer  edge 
of  the  boundary  layer,  respectively.  These  equations  fulfill  the  require¬ 
ments  of  no  slip  at  the  wall  as  well  as  prescribing  the  streamwise  distri¬ 
bution  of  wall  heat  flux.  The  outer  edge  velocity  ug  and  static  temperature 

T  are  obtained  from  the  inviscid  flow  calculation  of  Step  1  and  must  be 
e 

consistent  with  the  streamwise  distribution  of  static  pressure,  p(x). 


Before  a  solution  to  the  system  of  equations  defined  by  Equations 
5  through  13  is  possible,  the  form  of  the  turbulent  eddy  viscosity  must  be 
specified.  The  shear  stress  in  a  conventional  turbulent  wall  boundary  layer 
is  treated  herein  by  the  use  of  a  two-layer  inner-outer  model  based  on  the 
Prandtl  mixing-length  hypothesis.  That  is 


=  f2 


3U 


ay 


04) 


where  the  mixing  length  i  is  given  by  Escudier  (1965)  as 


n  =  *y  for  0  <  y  <  —  and 
i  =  A  6  for  —  <  y  <  6. 

K 


(15) 


Patankar  and  Spalding  (1968)  have  recommended  that  the  numerical  constants  be 
taken  as  k  =  0.435  and  A  =  0.09.  The  edge  of  the  boundary  layer  6  is  defined 
as  the  point  where  the  local  velocity  is  equal  to  ninety-nine  percent  of  that 
in  the  free  stream. 


In  the  near-wall  region,  where  both  laminar  and  turbulent  com¬ 
ponents  of  the  total  shear  stress  are  important,  Van  Driest's  (1956)  modifi¬ 
cation  to  the  mixing  length  has  been  employed.  The  effects  of  pressure 
gradient  and  heat  transfer,  which  are  commonplace  in  the  present  application, 
have  been  accounted  for  in  the  mixing-length  formula  by  Cebeci  and  Smith 
(1974). 
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This  basic  two-equation  model  for  conventional  turbulent  wall 
boundary  layers  has  been  extended  to  handle  the  case  of  wall  jets  in  the 
following  manner.  As  long  as  the  wall  jet  boundary  layer  exhibits  a  local 
jet  maximum  in  its  velocity  profile  the  remnant  of  the  main-stream  boundary 
layer  and  the  wall  jet  which  constitute  the  combined  layer  are  considered  to 
be  separate  entities.  Hence  the  two-equation  model  is  applicable  to  each 
region  independently.  The  eddy-viscosity  profile  is  made  continuous  between 
eddy- viscosity  maxima  by  fitting  a  cosine  fairing  between  these  maxima  as 
performed  by  Dvorak  (1973).  This  completes  the  set  of  closure  assumptions 
that  are  necessary  before  a  solution  is  possible. 

The  differential  equations  for  the  conservation  of  mass,  momentum 
and  energy  are  expressed  in  finite-difference  form  using  three-point  central 
differencing  in  the  y-direction  and  three-point  upstream  differencing  in  the 
x-direction.  In  this  manner  each  of  the  momentum  and  energy  equations  breaks 
down  into  a  system  of  linear  algebraic  equations  in  tridiagonally-banded  form 
that  is  solvable  by  rapid  efficient  means.  Once  the  inflow  boundary  con¬ 
ditions  are  specified  the  remainder  of  the  flow  field  is  solved  by  marching 
in  a  downstream  manner  from  the  near-slot  region  to  the  end  of  the  nozzle 
ring  element. 

2.3  Iteration 

Unfortunately,  as  in  the  case  of  many  coupled  two-step  solution 
procedures  for  aerodynamic  problems,  the  design  of  a  film-cooled  entraining 
diffuser  involves  iteration.  This  is  necessitated  by  the  fact  that  flow 
development  depends  on  factors  such  as  pressure  gradient  and  radiation  shape 
factors.  These  quantities  are,  in  part,  functions  of  geometry,  the  details 
of  which  may  not  be  known  until  the  flow  field  is  computed.  In  a  case  where 
a  given  degree  of  film-cooling  protection  is  desired,  for  example,  the  length 
of  nozzle  ring  element  that  is  acceptable  will  not  become  apparent  until  the 
computations  are  carried  out.  Only  then  will  it  be  seen  where  the  ring 
temperature  exceeds  some  prespecified  uppc*'  limit.  However,  before  such 
calculations  can  be  performed  the  pressure  gradient  that  is  imposed  on  the 
boundary  layer  must  be  defined.  This  depends  on  the  unknown  ring  length. 
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Consequently,  it  becomes  necessary  to  guess  at  this  unknown  length  and  to 
check  its  validity  a  posteriori. 

Other  iterations  may  be  necessary  in  the  event  that  the  sum  of  the 
computed  ring  lengths  exceeds  some  overall  limiting  length.  It  would  then  be 
necessary  either  to  alter  the  criterion  which  governs  ring  length  or  to 
adjust  the  slot  width  in  an  effort  to  accomplish  the  desired  end  result  in  a 
slightly  different  manner. 

A  flow  chart  which  illustrates  all  possible  iterations  is  shown  in 
Figure  5.  Depending  on  the  objectives  of  the  design,  the  constraints  and  the 
quality  of  first  guesses,  it  may  be  possible  to  eliminate  one  or  more  of  the 
loops  shown. 

3.0  SAMPLE  DESIGN  PROBLEM 

In  order  to  illustrate  the  design  procedure  a  component  of  equip¬ 
ment  recently  designed  at  DRES,  which  incorporates  a  film-cooled  entraining 
diffuser,  will  serve  as  an  example.  The  details  of  this  apparatus  a^e 
unimportant  as  far  as  this  exercise  is  concerned.  Suffice  it  to  say  that  the 
objectives  of  the  design  are 

(a)  to  achieve  a  wall  film-cooling  efficiency  in  excess 
of  95  percent,  and 

(b)  to  accomplish  significant  exhaust  stream  diffusion 
subject  to  the  following  constraints: 

i)  The  film-cooled  diffuser  must  be  compatible  with  a 
duct  of  28-inch  diameter.  This  duct  is  of  unspecified 
upstream  configuration  but  can  be  considered  a 
black  body  cooled  to  an  efficiency  of  95  percent, 
ii)  The  overall  length  of  the  structure  must  be  held  to 
36  inches  (measured  in  the  axial  direction  from  the 
first  slot  exit)  with  the  exit  diameter  not  greater 
than  33  inches  due  to  restrictions  on  available 
space. 
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FIGURE  5  :  FLOW  CHART  FOR  THE  DESIGN  OF  A 
MULTIPLY-SLOTTED  ASSEMBLY. 
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iii)  The  film-cooled  diffuser  is  to  be  constructed  of 
1/1 6- inch  mild  steel  whose  surface  is  treated  to 
exhibit  an  emissivity  of  unity, 

iv)  The  flow  of  hot  gas  is  turbine  exhaust  at  a  uniform 
velocity  of  151.9  fps,  a  temperature  of  1048°R  and 
an  effective  gray-body  emissivity  of  0.05. 

v)  Ambient  pressure  and  temperature  are  13.40  psia 
and  520°R,  respectively. 

Step  1 :  As  an  initial  guess  it  will  be  assumed  that  the  structure  con¬ 
sists  of  three  nozzle  ring  elements  which  overlap  so  as  to  create  three 
entraining  slots,  each  of  3/4-inch  width.  In  combination  with  the  ring 
material  thickness  this  gives  an  overall  outside  diameter  of  32-7/8  inches 
at  the  exit.  Note  that  this  dimension  is  just  under  the  maximum  allowable 
diameter  of  33  inches.  The  resulting  diffusion  will  therefore  approach 
the  maximum  possible  under  these  geometric  constraints. 

In  order  to  calculate  the  rate  of  entrainment  into  each  of  the 
three  slots  an  assumption  must  be  made  about  entrance  loss  coefficients. 

These  can  be  calculated  from  formulae  found  in  any  standard  text  on  elementary 
fluid  mechanics,  such  as  that  by  John  and  Haberman  (1971).  For  the  geometry 
of  the  present  application  the  loss  coefficients  for  slots  1,  2  and  3  have 
been  calculated  to  be  0.176,  0.150  and  0.150,  respectively.  These  co¬ 
efficients  are  typical  of  a  well-rounded  entrance  with  wiggle  stripping 
located  in  the  overlap  region  in  order  to  maintain  a  constant  slot  spacing. 

The  loss  coefficient  is  somewhat  higher  for  slot  number  one  due  to  a  15- 
degree  bend  in  the  overlap  region  (as  imposed  by  the  apparatus  upstream  of 
the  diffuser)  and  the  presence  of  fastening  clips  which  are  required  to 
secure  the  diffuser  to  the  upstream  duct. 

Results  of  the  entrainment  calculation  are  presented  in  Figure  6. 

As  shown,  the  mean  velocities  in  slots  1,  2  and  3  are  42.5,  35.5  and  25.8 
feet  per  second,  respectively. 

Step  2:  Having  computed  these  velocities  it  is  now  possible  to  calculate 
the  flow  development  downstream  of  each  slot.  An  initial  guess  must  be 
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FIGURE  6  :  RATES  OF  ENTRAINMENT  INTO  THE  THREE-RING 
ASSEMBLY  AS  CALCULATED  BY  PROGRAM 
'ENTRAIN' . 
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made  at  the  ring  lengths  j^,  i2  anc'  *3*  Experience  has  shown  that  the 
first  ring  is  usually  the  shortest  and  that  ring  length  increases  to  a 
maximum  somewhere  in  the  middle  of  the  structure  and  then  decreases  toward 
the  downstream  end.  This  trend  is  the  outcome  of  several  factors: 

i)  Entrainment  velocity  decreases  in  the  downstream  direction 
due  to  a  decelerating  main  stream  and  an  increasing  area 
ratio. 

ii)  Radiative  heat  transfer  loads  decrease  as  the  distance  from 
the  exhaust  duct  increases. 

iii)  The  build-up  of  cooling  layers  tends  to  improve  cooling 
performance  at  downstream  locations. 

As  a  first  guess  the  ring  lengths  i2  and  will  be  10,  14  and  12 

inches,  respectively,  to  give  an  overall  length  L  of  36  inches  as  given 
in  the  list  of  constraints.  The  pressure  gradient  along  each  ring  can  now 
be  calculated  using  the  ring  lengths  above  and  static  pressures  in  Figure 
6.  As  well,  now  that  the  geometry  has  been  specified,  the  distribution  of 
radiation  shape  factor  can  be  computed.  A  fairly  simple  approach  to  this 
calculation  is  outlined  in  Appendix  B  where  a  computer  program  'SHAPE'  is 
described.  The  results  of  'SHAPE'  for  this  example  are  summarized  in 
Figure  7.  Two  shape  factors  are  listed  for  each  position  along  the  ring 
structure  at  1-inch  axial  intervals.  One  shape  factor,  F^_q,  relates  to 
radiative  heat  transfer  from  the  upstream  duct  to  a  ring  element  of 
infinitesimal  width  at  the  axial  station  in  question.  The  other  shape 
factor,  Fcb,  relates  to  radiative  heat  transfer  from  this  ring  element  to 
the  atmosphere  as  viewed  through  the  exit  end  of  the  diffuser  assembly. 

Before  a  boundary-layer  calculation  is  possible  assumptions 
must  be  made  about  the  velocity  profile  in  each  of  the  slots  and  in  the 
exhaust  duct  just  upstream  of  the  first  slot.  These  assumptions  must  be 
made  by  the  designer  based  either  on  his  knowledge  of  the  upstream  flow 
or  on  experimental  data.  With  the  aid  of  the  "law  of  the  wall"  and 
Coles'  (1956)  "law  of  the  wake"  realistic  velocity  profiles  can  be 
constructed  once  basic  information  about  the  boundary  layers  is  provided. 
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FIGURE 


RADIATION  SHAPE  FACTORS  FOR  A  TAPERED  DUCT 
*ITH  CAPPED  t NOS  A  AND  B 
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7  :  RADIATION  SHAPE  FACTORS  FOR  THE  THREE-RING 
ASSEMBLY  AS  CALCULATED  BY  PROGRAM  'SHAPE'. 
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In  addition  to  satisfying  the  above  laws  flow  in  the  slot  must  meet  two 
other  criteria.  Firstly,  since  the  slot  flow  is  one  of  low  Reynolds 
number.  Coles'  profile  parameter  n  is  a  function  of  Reynolds  number 
(based  on  the  momentum  deficit  thickness).  Secondly,  the  velocity 
profile  must  be  constructed  such  that,  when  integrated  across  the  slot, 
it  yields  the  same  mean  velocity  as  calculated  in  the  entrainment  analysis 
of  Step  1.  These  constraints  necessitate  iterative  procedures  in  setting 
up  the  slot  velocity  profiles.  Parameters  used  to  create  the  velocity 
profiles  in  the  present  example  are  summarized  in  Table  I  below. 


Flow 

Maximum  Streamwise 

Velocity  u„ 

J  e 

fps 

Boundary-Layer 
Thickness  6 

inches 

Coles'  Profile 
Parameter  n 

main  stream 

151.9 

2.0 

0.550 

slot  1 

48.5 

0.25 

0.141 

slot  2 

40.3 

0.25 

0.080 

slot  3 

29.3 

0.25 

0.000 

TABLE  I:  Boundary-Layer  Parameters  used  to  Define  the 
Main  Stream  and  Slot  Flows 

Results  of  the  boundary-layer  calculation  are  shown  in  Figures  8 
and  9.  All  distances,  velocities  and  temperatures  are  in  units  of  inches, 
fps  and  °R,  respectively.  The  vertical  scale  on  each  plot  has  been  exag¬ 
gerated  by  a  factor  of  two  in  relation  to  the  horizontal  scale.  The  finite 
grid  used  to  represent  the  flow  field  consists  of  36  equal-spaced  axial 
stations  and  approximately  100  variable-spaced  radial  stations.  About  eight 
minutes  of  central  processing  time  were  required  on  an  IBM  360  time-sharing 
system.  Double  precision  (eight-byte  variables)  was  required. 

Figure  8  shows  velocity  profiles  at  every  second  axial  station. 

For  the  low  entrained  air  velocities  which  are  typical  of  the  present 
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FIGURE  8  :  VELOCITY  PROFILES  IN  THE  THREE-RING 
ASSEMBLY  AS  CALCULATED  BY  PROGRAM 
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study  the  jets  are  seen  to  degenerate  rather  rapidly.  This  is  more  pronounced 
at  downstream  slots  where  the  velocity  of  the  entrained  stream  becomes  pro¬ 
gressively  smaller. 

Temperature  profiles  and  isotherms  are  shown  in  Figure  9.  The 
temperature  profiles  show  clearly  the  degeneration  of  the  zero  temperature- 
gradient  zone  in  the  near-slot  region.  One  important  feature  of  the  plot  that 
is  not  evident  due  to  the  scale  chosen  is  that  temperature  profiles  near  the 
slot  exhibit  a  large  gradient  immediately  adjacent  to  the  wall.  This  is  so 
because  radiative  heat  transfer  causes  the  wall  to  be  at  a  higher  temperature 
than  the  neighbouring  fluid  above  it.  This  feature  is  more  apparent  on  the 
isotherm  plot  in  that  the  isotherm  adjacent  to  the  wall  exhibits  a  discon¬ 
tinuity  and  change  in  direction. 

Another  observation  of  interest  is  that  isotherms  are  mildly  discon¬ 
tinuous  in  the  near-slot  region.  These  discontinuities  occur  because  the 
eddy-viscosity  model  (responsible  for  describing  the  turbulent  mixing  process) 
changes  once  the  local  wall  jet  maximum  in  the  velocity  profile  disappears. 

4.0  EXPERIMENTAL  VERIFICATION 

In  order  to  validate  the  design  procedure  presented  in  2.0,  it  was 
decided  to  build  and  test  the  structure  described  in  the  example  of  3.0.  A 
detailed  assembly  drawing,  shown  in  Figure  10,  was  prepared  based  on  computed 
slot  widths  and  ring  lengths.  A  photograph  of  the  manufactured  structure 
appears  in  Figure  11.  It  is  seen  to  consist  of  three  overlapping  mild-steel 
rings  fastened  together  with  steel  wiggle-type  stripping  in  order  to  maintain 
a  3/4-inch  slot  spacing.  Metal  tubing  of  circular  cross  section  was  machined 
and  welded  to  the  leading  end  of  each  ring  so  as  to  create  smooth  and  well- 
rounded  slot  entries.  The  inner  surface  of  each  ring  element  was  treated  to 
achieve  the  highly-absorbing  finish  desired  for  the  present  application. 

Although  it  is  recognized  that  detailed  experimental  validation 
should  include  hot-wire  anemometry  and  resistance- thermometry  data  in  order 
to  verify  that  velocity  and  temperature  fields  are  as  predicted,  a  much  simpler 
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FIGURE  10  :  ASSEMBLY  DRAWING  OF  THE  THREE-RING  ASSEMBLY. 
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approach  will  be  taken  here  due  to  limited  time  and  resources.  For  the 
present,  it  is  proposed  that  the  design  procedure  be  checked  by  two  sets  of 
simple  measurements.  In  the  first  set  stagnation  probes  and  static  pressure 
taps  will  be  used  to  evaluate  entrance  loss  coefficients  and  entrained  air 
velocities.  The  second  set  of  measurements  will  use  thermocouples  to  measure 
skin  temperature  at  various  streamwise  locations  along  the  structure.  If 
satisfactory  results  are  not  obtained  it  may  be  necessary  to  revert  to  more 
elaborate  means  of  instrumentation. 

It  is  intended  to  perform  the  experiments  in  the  DRES  film-cooling 
facility  shown  schematically  in  Figure  12.  The  facility  consists  of  an  Orenda 
gas  turbine  engine  which  supplies  hot  gas  to  a  test  section  downstream  of  a 
flow  control  bypass  vane. 

5.0  CONCLUSIONS 

An  iterative  design  procedure  for  film-cooled  entraining  diffusers  has 
been  developed  and  results  have  been  presented  for  a  specific  application.  A 
test  article  based  on  this  design  has  been  fabricated  and  experiments  with  this 
apparatus  will  be  conducted  in  the  DRES  film  cooling  facility.  Results  of 
these  tests  will  be  reported  separately  at  a  later  date. 
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FIGURE  12  :  THE  ORES  FILM-COOLING  FACILITY. 
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APPENDIX  A:  PROGRAM  'ENTRAIN' 


The  system  of  equations  defined  by  Equations  1  through  5  can  be 
rearranged  so  that  the  velocity,  static  pressure  and  temperature  at  the 
exit  plane  of  any  nozzle  ring  element  are  given  by 
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The  above  method  is  used  in  the  iterative  solution  method 
employed  in  program 'ENTRAIN'.  A  documented  listing  of  the  program  appears 
in  Figure  A-l. 
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FIGURE  A-l :  A  DOCUMENTED  LISTING  OF  PROGRAM  'ENTRAIN' 
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APPENDIX  B:  PROGRAM  'SHAPE1 

Two  modes  of  radiative  heat  transfer  are  present  in  the  diffuser: 
radiation  between  the  diffuser  walls  and  sources  both  upstream  and 
downstream  of  the  diffuser  assembly,  and  radiation  from  the  hot  exhaust 
gas  to  the  walls  of  the  diffuser  assembly.  The  combined  effect  of  these 
modes  of  heat  transfer  is  considered  to  be  the  simple  additive  result. 

This  approach  is  justified  since  the  gas  is  largely  transparent. 

Consider  the  first  mode  of  radiation  described  above.  The 
diffuser  is  of  relatively  complex  geometry  with  regard  to  the  calculation 
of  radiation  shape  factors.  However,  since  one  of  the  goals  of  the 
present  application  is  to  ensure  a  relatively  constant  wall  temperature 
via  the  use  of  air  film  cooling,  radiative  heat  transfer  between  neighbouring 
regions  within  the  diffuser  is  small.  Clearly,  the  significant  exchange 
of  radiant  energy  is  between  the  walls  of  the  diffuser  and  both  the 
upstream  duct  and  downstream  ambient  surroundings.  With  this  in  mind 
the  calculation  of  shape  factors  is  simplified  considerably.  Although 
the  diffuser  consists  of  a  series  of  staged  cylindrical  rings,  for  the 
purpose  of  computing  these  factors,  it  will  be  approximated  by  a  conical 
surface  whose  cross-sectional  radius  varies  linearly  in  the  axial  direction 
from  that  of  the  exhaust  duct  to  that  of  the  largest  cylindrical  ring  as 
illustrated  in  Figure  B-l.  The  most  significant  errors  incurred  by  this 
approximation  are  for  the  near-slot  region  where  the  step  created  by  the 
slot  tends  to  shield  the  wall  to  a  higher  degree  than  this  approach 
would  indicate.  Since  the  amount  of  taper  is  small  these  errors  are 
considered  to  be  negligible. 

From  the  point  of  view  of  radiation  the  upstream  duct  is  well 
represented  by  a  disk  of  the  same  radius,  temperature  and  emissivity 
situated  at  the  entrance  of  the  diffuser.  Likewise,  a  perfectly  absorbing 
disk  at  ambient  temperature  located  at  the  exit  plane  simulates  the  cool 
ambient  surrounding.  These  disks  are  labelled  A  and  B  respectively  in 
Figure  B-l. 

The  shape  factor  from  disk  A  to  an  arbitrary  ring  element  C  of 
length  2dx  is  simply  the  difference  between  the  shape  factor  from  disks  A  to 
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FIGURE  B-l :  SIMPLIFIED  DIFFUSER  FOR  THE  PURPOSE  OF 
CALCULATING  RADIATION  SHAPE  FACTORS, 
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Cj,  and  that  from  disks  A  to  C2.  According  to  McAdams  (1942)  the  shape 
factor  from  disk  A  to  disk  C.,  centred  on  the  same  axis  a  distance  x  apart,  is 


A-C, 


1 

2a2 


+  c.2  -  V(x2  +  a2  +  c.2F-  4a2ci<2j 


where  a  and  c-  are  the  radii  of  disks  A  and  C.. ,  respectively.  The  heat 
transfer  per  unit  area  of  the  conical  ring  element  described  above  is  then 


AA 

^A-C  =  oeAEC  ^A4  ”  V>  ^a-C 

where  FA-C  =  ^Tci+^rdT  *  ~T 
-V((x  -  dx)2  +  a2  +  Cl2)2  -  4a2c, 


[(x  -  dx)2  +  a2  +  C] 


2]  -  [(x  +  dx)2  +  a2 


2 


+  C2 


2 


-V((x  +  dx)2  +  a2  +  c22)2  -  4a2c22]J  . 

Here  e  is  gray-body  emissivity  and  T  is  absolute  temperature.  Taking  the 
limit  of  the  above  expression  as  dx  tends  to  zero,  with  the  aid  of  L'Hopital's 
rule,  yields 


qA-C  "  aeAEC^TA4  "  V> 


x2  +  a2  +  c2 
+  a2  +  c2)2  -  4a2c2 ' 


where  a  is  Boltzmann's  constant. 


A  similar  analysis  for  the  radiative  heat  transfer  from  the  ring 
element  to  disk  B  gives 


qC-B  =  aeBeC  (V  ‘  V)  ^o" 


2c 


(L  -  x)2  +  b2  +  c2 


V((L  -  x)2  +  b2  +  c2)2  -  4b2c2' 


-  1 


The  second  mode  of  radiative  heat  transfer,  namely  that  from  the  hot 
exhaust  gas  to  the  diffuser  walls,  is  calculated  by  selecting  a  suitable 
effective  gray-body  emissivity  for  the  hot  gas  and  by  assuming  a  shape  factor 
near  unity  since  the  column  of  gas  almost  completely  fills  the  interior  of  the 
diffuser.  Effective  gray-body  emissivity  for  industrial  and  power  plant 
exhaust  gases  having  various  concentrations  of  water  vapour,  carbon  dioxide 
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and  carbon  particulates  is  given  by  Holman  (1976). 

A  simple  computer  program  'SHAPE'  which  calculates  the  area- 
normalized  shape  factors  above  appears  in  Figure  B-2. 
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c  *  *••****«*•*•**  ***************  SHAPE. 

SHAPE 

CO»M*****OMM*«M*****»*«**  SHAPfc 

c 

c 


C  The  purpose  of  THIS  program  is  to  calculate  radiation  shape 

C  FACTORS  BETWEEN  VARIOUS  SURFACES  IN  A  CONICAL  DUC T  WITH  CAPPED 

L  ENDS,  in  particular  the  fulluwing  shape  factors  are  desired# 

c 

C  1)  THAT  HETnEFN  the  SMALL  CAP#  A.  AND  A  RING  ELEMtNl,  C#  UF 

C  INFINITESIMAL  WIDTH#  ANO 

C 

C  2)  THAT  ttE  T  «F  E  n  THE  RING  ELEMENT  DESCRIBED  AHOVE  AND  THE 

C  LARGE  LAP#  H. 

c 

C  A  IS  The  RADIUS  UF  THE  SMALL  CAP  (INCHES). 

C 

L  B  IS  THt  RADIUS  OF  THE  LARGE  CAP  (INCHES). 

C 

C  XL  IS  THE  AXIAL  LtNGTH  UF  THE  CONICAL  DUCT  (INCHES),  ' 

C 

C  DX  IS  THE  AXIAL  DISTANCE  BETWEEN  POINTS  FOR  WHICH  THE  SHAPE 

C  FACTOR  IS  DESIRED  (INCHES). 

C 

C 

t  INITIALIZATION  AND  DATA  INPUT. 

C 

I U  T  =h 
tN  =  5 

RE  AD ( I N# 10)A, B#XL.DX 
10  FORMAT («F 10. S) 

KOUNTsl . 00001 *( XL/DX) 

wR I TE ( I UT  #  20 ) 

20  FORMAT (//////«1X# ’RAD1AJ ION  SHAPE  FACTORS  FOR  A  TAPERED  OUCT'// 

1  SOX#  *  W I TH  CAPPED  ENOS  A  AND  B  * / /// 

2  USX, 'AXIAL  RADIUS  SHAPE  SHAPE*/ 

3  43X, fPUSI T ION* # lax, 'FACTOR  F AC T OR ' //66X # • F  F' 

u  /aax, 'Inches  Inches  a-c  C-b'//) 

c 
c 

c  THt  CALCULATIONS  FOR  SHAPE  FACTOR  BEGIN. 

C 

DO  90  J* 1 # KOUNl 
IF (KOUNt-J) JO. 30. UO 
3U  X=XL-0.001 *DX 
GO  TO  50 
ao  XsJ*DX 
SO  CSA+X/XL*(B-A) 

C 

L 

C  THE  SHAPE  FACTOR  FRUM  A  TO  C  IS  COMPUTED  BELOW. 

C 

PART 1=X«XFA*A*C«^ 

PAKT?SSUR  T (PART  1  APART  1-9. «A*A*C*C  ) 

SHPAC=x/2./C* (PARTI /PART 2-1 .01 
C 
L 

C  THE  SHAPE  FACTOR  FROM  C  TU  H  IS  COMPUTED  BELOW. 

C 

PART  5=(XL-X  )  *  UL-X>^tJ«<B*CAt 

PART «»SURT (PART  3 APART J-a. «B*B*C*C ) 

SMPCHs ( XL-X )/2./C*( PART  3/PART  a- 1 .0) 

C 

C 

C  RESULTS  ARt  PRINTED  BELOw. 

C 

IF (KOUNT-J)60#60, 70 
bO  XrXL 

70  wRI TE ( IOT,BO)X,C#SHPAC#SHPCB 
80  FORMAT (  u  1  x  #  uF 10. S) 

90  CONTINUE 
C 
L 

STOP 
E  NO 


FIGURE  B-2:  A  DOCUMENTED  LISTING  OF  PROGRAM  'SHAPE'. 
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